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Beyond the quantum Markov approximation and the weak coupling limit, we present a general 
theory to calculate the geometric phase for open systems with and without conserved energy. As 
an example, the geometric phase for a two-level system coupling both dephasingly and dissipatively 
to its environment is calculated. Comparison with the results from quantum trajectory analysis is 
presented and discussed. 
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I. INTRODUCTION 

Since Berry's discovery that a state of a quantum 
system can acquire a phase of purely geometric origin 
when the Hamiltonian of the system undergoes a cyclic 
and adiabatic change, there have been numerous pro- 
posals for generalizations, including the geometric phase 
for nonadiabatic, noncyclic, and nonunitary evolution 
the geometric phase for mixed states 0, Q , the geometric 
phase in systems with quantum field driving and vacuum 
induced effects ^ , as well as the geometric phase in cou- 
pled bipartite systems^]. 

Recently, the attention has turned to studying geo- 
metric phase in open systems, this is motivated by the 
fact that all realistic system are coupled, at least weakly, 
to their environment. From the perspective of applica- 
tion, the use of geometric phases in the implementation 
of fault-tolerant quantum gates 0, H, 0, requires the 
study of geometric phases in more realistic systems. For 
instance, the system which carries information may de- 
coher from a quantum superposition into statistical mix- 
tures, and this effect, called decoherence, is the most im- 
portant limiting factor for quantum computing. 

The study of geometric phase in open systems may be 
traced back to 1980s, when Garrison and Wright ^^^^ 
touched on this issue by describing open system evolu- 
tion in terms of non-Hermitian Hamiltonian. This is a 
pure state analysis, so it did not address the problem of 
geometric phases for mixed states. Toward the geometric 
phase for mixed states in open systems, the approaches 
used include solving the master equation of the system 
[l^lT^ l ^ f islll^ . employing a quantum trajectory anal- 
vsis [ivHlSl or Krauss operators the perturba- 

tive expansions |^. Some interesting results were 
achieved that may be briefly summarized as follows: non- 
hermitian Hamiltonian lead to a modification of Berry's 
phase 10,1131, stochastically evolving magnetic fields pro- 
duce both energy shift and broadening |2l|. phenomeno- 
logical weakly dissipative Liouvillians alter Berry's phase 
by introducing an imaginary correction |l4j or lead to 
damping and mixing of the density matrix elements |l5j| . 



However, almost all these studies are performed for dis- 
sipative systems under various approximations, thus the 
representations are approximately applicable for systems 
whose en ergy is not conserved. Quantum trajectory anal- 
ysis ^3, Il8l | that uses the quantum jump approach is 
available for open systems with conserved energy. Its 
starting point, however, is the master equation, a result 
within the quantum Markov approximation and in the 
weak coupling limit. Beyond the quantum Markov ap- 
proximation and the weak coupling limit, the geometric 
phase of a two-level system with quantum field driving 
was analyzed |12 , where the whole system (the two- level 
system plus the quantum field) was assumed subject to 
dephasing, the subsystem(two-level system in that paper) 
may still acquire geometric phase even when the whole 
system in its pointer states. This is an ideal situation to 
show the vacuum effects on the geometric phase of the 
subsystem, as well as the decoherence effects on the geo- 
metric phases regardless of its feasibility of experimental 
realization. However, beyond the Markov approximation 
and the weak coupling limit, the geometric phase for a 
dissipative system remains untouched. In this paper, we 
will deal with the geometric phase in open systems, be- 
yond the Markov approximation and weak coupling limit. 

The structure of this paper is organized as follows. In 
Sec. II the exact solution and calculation of the geometric 
phase of a system dephasingly coupled to its environment 
is presented, an example to detail the representation and 
a discussion on physical realization are given in Sec. HI. 
In Sec. IV, we present an example to show the calculation 
of geometric phases in dissipative systems. Finally we 
conclude in Sec.V. 



II. GEOMETRIC PHASE IN DEPHASING 
SYSTEMS: GENERAL FORMULATION 

In this section, we investigate the behavior of geomet- 
ric phase of a quantum system subject to decoherence. 
In order to make a comparison with the results from the 
quantum jump approach, we consider the quantum sys- 
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tern without any field driving except the environment. 
So, it is not directly relevant to our previous study [l^. 
The environment that leads to decoherence may origi- 
nate form the vacuum fluctuations or the background 
radiations. We restrict ourselves here to consider the 
case where the system-environment coupling Hi and the 
free system Hamiltonian Hs commute. This is the sit- 
uation of dephasing and the exact analytical dynamics 
may be obtained. On the other hand, the evolution of 
a system with such properties may be described by the 
master equation when the Markov approximation and the 
weak coupling assumption apply, this kind of decoherence 
would not change the geometric phase of the quantum 
system by the quantum jump approach [T^ IT^ . However, 
as you will see, this is not the case from the viewpoint of 
interferometry considered in this paper. 

We consider a situation described by a Hamiltonian of 
the form 

H^Hs + Hb+Hj, (1) 

where Hs describes the free Hamiltonian of the system, 
Hb stands for the Hamiltonian of the environment, and 
Hi represents the system-environment couplings. The 
environment and the system Hamiltonian may be taken 
arbitrary but with constraints [Hs, Hi] — 0. Let us sup- 
pose that the interaction Hamiltonian Hi has the form 
(setting h — 1) 

Hi = Y,Xrn{rl + r^), (2) 

m 

where the Xjn,{m = 1,2,...,M) are the system oper- 
ators satisfying [Hs,Xjn] = 0, and the represent 
environment operators that may take any form in gen- 
eral. Commutation relation [Hs,Xra] = enables us 
to write the time evolution operator for the whole sys- 
tems(system-l-environment) as 

Uit) = e-'^* = e-*^-* UnME„.){Ern\ (3) 

m 

with Um{t), a function of environment operators satisfy- 
I 



Here, $3 is the geometric phase for the system under- 
going non- unitary evolution |23|. T denotes a time after 
that the system completes a cyclic evolution when it is 
isolated from any environments. Taking the effect of en- 
vironment into account, the system no longer undergoes 
a cyclic evolution, but we still keep the T as the period 
for calculating the geometric phase in this paper. The 
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He,m = HB + J2ei:{T^ + ri). (4) 

n 

Here, \E^) stands for the eigenstate of Hs with eigen- 
value i?m|2a|j while e™ denotes the eigenvalue of X„ cor- 
responding to eigenstate For a specific F^, Um{t) 
may be expressed in a factorized form, which will be 
shown later through the spin-boson model. Furthermore, 
we assume that the environment and the system are ini- 
tially independent, such that the total density operator 
factorizes into a direct product, 

piO) ^ ps{0) ® Pb{0) = 

J2pmnmE,n){En\<»PB{0). (5) 

mn 

At time t, the reduced density operator of the system is 
given by 

Ps{t) = TrBiU{t)ps{0)(E>PBiO)U\t)) 

= ^Pmn(0)e-*(^"-^")*|i?™)(i?„|F™„(t), (6) 

mn 

where F^n{t) is defined as TrB[Um{t)pB{0)Ul{t)]. Eq.© 
shows that the diagonal elements of the reduced density 
matrix Pmm are time-independent, while the off-diagonal 
elements evolve with time involving contributions from 
the environment-system couplings, it at most cases would 
lead to decaying in the off-diagonal elements, and even- 
tually results in vanishing of these matrix elements. Now 
we turn to study the geometric phase of the open sys- 
tem. For open systems, the states in general are not 
pure and the evolution of the system is not unitary. For 
non-unitary evolutions, the geometric phase can be cal- 
culated as follows. First, solve the eigenvalue problem for 
the reduced density matrix p{t) and obtain its eigenvalues 
ek{t) as well as the corresponding eigenvectors \ipk{t)); 
secondly, substitute ek{t) and into 



(7) 

I 

geometric phase Eq. Q is gauge invariant and can be 
reduced to the well-known results in the unitary evolu- 
tion, thus it is experimentally testable. The geometric 
phase factor defined by Eq.© may be understood as a 
weighted sum over the phase factors pertaining to the 
eigenstates of the reduced density matrix, thus the detail 
of analytical expression for the geometric phase would de- 



= Re{arg(^ V^;(0)^(^,(0)|^..(r))e-r<'^'=(*)l^/^*l'^'=(*»'^*)}. 
fc 
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pend on the digitalization of the reduced density matrix 
Eq.®. 



III. GEOMETRIC PHASE IN DEPHASING 
SYSTEM: EXAMPLE 

To be specific, we now present a detailed model to illus- 
trate the idea in Sec. 11. The system under consideration 
consists of a two-level system coupling to its environment 
with coupling constants {gi}, the Hamiltonian which gov- 
erns the evolution of such a system may be expressed as 



H = 



LO 



(|e)(e| ~ \g){g\) + \{\e){e\ - \g) {g\)Y. g^{a\ + a.) 



E 



where ai are the creation and annihilation opera- 
tors of the environment bosons, and |e), \g) denote the 
excited and ground states of the two-level system with 
Rabi frequency uj. This Hamiltonian corresponds to 
Xm = X = \{\e){e\ - \g){g\), and r„ = T = J^idi^i in 
the general model Eq.|21). Generally speaking, the choice 
of the coupling between the system and the environment 
determines the effect of the environment. For exam- 
ple, the choice of the system operator Xm that do not 
change the good quantum number of Hs when they op- 
erate on the eigenstates of Hs would result in dephasing 
of the system, but not energy relaxations. The system- 
environment coupling taken in this section is exactly of 
this kind. 

By the procedure presented above, the reduced density 
matrix in basis {|e), \g)} for the open system follows [2J|, 



PS 



cos 



2 e 



isin6lFi2(i) 



isine'i^2i(i) sin^ 



(9) 



where an initial state of (cos||e) + sin||g)) (g) \0)b for 
the total system was assumed in the calculation, and 



Fuit) = F2i{t) = F(t) = 



(10) 



with rij{t) — 4|^p(l — cosujjt), and \0)b denoting the 
vacuum state of the environment. Some remarks on the 
reduced density matrix are now in order. For any j, 
rij{t) > 0, so as t tends to infinity (with respect to the 
system's coherence time), F{t) tends to zero, this indi- 
cates that the off-diagonal elements would vanish on a 
long time scale with respect to the decoherence time, 
and hence the open system would not acquire geomet- 
ric phase when time is longer than the decoherence time. 
This is different form the results concluded in the pre- 
vious work, where the subsystem may acquire geometric 
phase even for the whole system in its pointer states . 
To calculate the geometric phase pertaining to Eq.®, 
we first write down the eigenstate and its corresponding 



eigenvalue for the reduced density matrix ps as, 
\e±{t)) = cose±(t)|e) -|-sine±(t)|.9). 



e±{t) - i(l±y^cos2 + sin2 0|^^(i)P), (11) 



with 
cos8±(i) 



sm0F{t) 



sin2e'|F(i)|2-|-4(e±(t)-cos2f)2 



(12) 



Clearly, for a closed system, namely gi = 0, F{t) = e 
the eigenvalue and corresponding eigenstates reduce to 



£±(^) — IjO, and cos9+(t) — cos ■ 



*, sine+(i) = 



sinf, cose_(t) = -sinf, smO-it) = cosfe*"*. Ah 
these together yield the well-known geometric phase 

(8) $3"^ = (1 + cos 6*). Eq.ini) and Eq.^ are the exact re- 
sults for the open two-level system, the geometric phase 
would depend on how F(t) varies with time. For a 
continuous spectrum of environmental modes with con- 
stant spectral density a{iu) = e, F{t) = e~"^*e~'>'* with 
7 = 27re|(7|2, where gi = g was assumed. Up to the first 
order in 7, the geometric phase at time T = 2tt/uj is given 
by 



$(1) 
g 



7r(l + cos 9) - — TT^ sin^ I 

UJ 



(13) 



This result can be easily understood as follows. The 
geometric phase factor for mixed states is defined as a 
weighted sum over the phase factors pertaining to the 
eigenstates of the reduced density matrix, the dephas- 
ing that leads to decaying in the off-diagonal elements 
would change the phase factors acquired by each eigen- 
state of the reduced density matrix, thus it modifies the 
geometric phase. This is different from the definition in 
the quantum jump approach 18], in which the problem 
of defining Berry's phase for mixed states was avoided 
by approaching the dynamics of open system from a se- 
quence of pure states, this leads to the result that the 
geometric phase is unaffected by dephasing, but it low- 
ers the observed visibility in any interference measure- 
ments. From the aspect of mixed state, the evolution 
of the system is among several trajectories with corre- 
sponding probabilities, so the geometric phase is defined 
as a weighted sum over the trajectories that the system 
undergoes. 

Now we are in a position to discuss the geometric phase 
acquired at time T — 2t:/lu by the two- level system. Sub- 
stituting Eqsimj and in to Eq.Q, we obtain 



cofi^ Q+{t)dt. 



(14) 



It is worthy to notice that Eq. (|14|l is the geometric phase 
beyond the quantum Markov approximation and the 
weak coupling limit. In this sense it provides us more 
insight into the geometric phase for dephasing system. 
The numerical results for Eq. H14|) are presented in fig- 
ure 1, where the dependence of the geometric phase is 
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FIG. 1: An illustration of the geometric phases of the two- 
level system coupling to the environment at time T — 2n/ui. 



illustrated as a function of the azimuthal angle 9 and 
the damping rate 7. The spectrum of environmental 
modes was taken to be cr(w) = e in this plot. Clearly, 
the two-level system acquires zero geometric phase with 
7 tends to infinity, this indicates that the two-level sys- 
tem acquires no geometric phase after the decoherence 
time. The representation in this paper may be used to 
analyze and estimate the error in the holonomic quan- 
tum computation due to decoherence [2^, in which the 
key error occurs within the degenerate subspace. The 
Hamiltonian that describes such a system reads H = 
J2^,j,a9?,J{ai + aa)\i){j\ + J2i'^iO'l"'i^ where the degen- 
erate energy was assumed to be zero, {\i),i — 1,...,N} 
denote the degenerate levels coupled to the environment 
with coupling constants {gi'j}- The Hamiltonian can be 
rewritten as H = I]^,„g^(aJ, -I- aa)\l3){f3\ + Y.^^iO'i"■^^ 

with an appropriate choice of \(3) — X^i^i This ex- 
actly is the case discussed in sections II and III. 

The open system effect in geometric phase may be ob- 
served with a combination of the engineering reservoir 
technique |27l and the Mach-Zchnder atom interferom- 
eter [23, [23, in which each of the arms consists of an 
atom in a dark state. The dark state can be realized in 
the atom-light system that consists of cesium atoms in- 
teracting with light resonant with the F = 3 ^ F =3 
transitions of the Dl 6^si/2 — > 6^Pi/2 line. It makes the 
dynamical phase negligible with respect to the geometric 
phase. A dephasing engineering reservoir in one arm of 
the interferometer, which may be simulated by variations 
of the light fields, would give a relative phase to the atom 
passing though the arm, the output interference pattern 
yields the geometric phase of the atom system. 

One of the key assumptions of our representation is the 
dephasing condition, i.e., [Hs,Hi] — 0. Using ground 
states as the qubits can make the formulation exact. 
Suppose now there is an small additional term in Hi, 
Hj = YmiYm + Y\^), that breaks the dephasing con- 
dition. Simple algebra shows that the transition proba- 
bihty between \Em) and |i?„) due to coupling YpiTp~\-T^p) 



is proportional to |7p(-E'„|'Kp|£'m)P/|-E™-£'„p, where 7p 
denotes the maximum of average values of (Fp -I- Fp. In 
the case of \^p{En\Yp\E„^)\ « \E,n ~ En\, the open sys- 
tem may be treated as a dephasing system, because the 
transition between any different eigenstates of the system 
may be ignored. The case where this transition could not 
be ignored will be discussed in the next section. 



IV. GEOMETRIC PHASE IN DISSIPATIVE 
SYSTEMS: EXACTLY SOLVABLE MODEL 



In this section, we will consider a spin-^ particle in- 
teracting with an environment formed by A'^ independent 
spins through the Hamiltonian 



1 ^ 



(fc) 



(15) 



fc=l 



where a^'^'' and ai, i — x,y,z denote Pauli operators 
for the environment and spin-i particle, respectively. 
gk,{k = l,2,...,iV) are coupling constants, term with 
A stands for the self-Hamiltonian of the particle. This 
model is interesting because the pointer states do not 
coincide with the eigenstates of the interaction Hamilto- 
nian, but can range from coherent sates to eigenstates of 
the system's Hamiltonian determining by the interplay 
between the self-Hamiltonian and the interaction with 
the environment. We will calculate the geometric phase 
gained by the particle beyond the Markov approxima- 
tion and the weak coupling limit. The dynamics govern 
by Hamiltonian Eg. ljlSII can be solved exactly by a stan- 
dard procedure's^, it yields the reduced density matrix 
of the particle as 



p(t) = (/ + p(t)-a)/2. 



(16) 



where p{t) is the polarization vector given by p{t) = 
f p{t, B)Tj{B)dB with n% ^ A"^ + B'^ and 



N 



, , A^ + B^ coa(2nBt) ,,B 
P.{t,B) = p,(0) ^ ' -py{0)—s\n{2^lBt) 



+ Pz{Q)-^SlT?{^Bt), 

Py{t, B) = py{Q) cos{2nBt) + [p^(o)i? - ApM] , 

Pz{t,B) = p^{0) -2 

"b 

2AB A 
+ pA0)-^sm^i^Bt)+Pyi0)—sm{2nBt). (17) 

To get this result, it is only required that the couplings gk 
of Eq. H15|) are sufficiently concentrated near their average 
value so that their standard deviation exists and is finite. 
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By rewriting the reduced density matrix p{t) in the we get the geometric phase[22 of the particle acquired 
form at time r, 

Pit) = Ai (0 1 V-i (t)} (V-i (i) I + A2 (t) I (t)) (V^2 (i) I , (18) 



(19) 



After simple manipulations, we arrive at 



,(^)(cos^cos^e-W-)-*(")) 



. e{0) . 0(t) 

sm sm — 

2 2 



2 2 
sin^^^)e''^ 



+ A2(r)(cos^cos^e'(*(-)-*(°)) 



• f(0) . 6'(t), _,ci 
sin^sm^)e '^]. 



(20) 



Here, 



Ai,2 = -V(l±|p(0)|)(l±b-(t)|), 



COS0{t) 



Pz(t) 



tan (/)(t) 



Py(^) 
Px{ty 



C 



dt 



-dt 



cos 9——dt 
dt 



The dependence of the geometric phase on the vari- 
ance sjv and system free energy A is complicated, we 
discuss here on two limiting cases stv S> A and sat <c A 
with a specific initial state Px{Q) = l,Py(0) = Pz(0) = 0. 
In sjv 3> A limit, the dynamics of the spin-i particle is so 

slow that its behavior should approach Px{t) = e"^* 
and Py{t) = Pz{t) — 0, which yields ^'g = because 
^ = in this limit with the initial state. In s^r <^ A 
limit, Eq.lEI) follows that. 



Px{t) 
Py{t) 



P.(0) 



Pz{t)^0, 



A 



\/2, 



SN 



(2 At 



3^T ^ 
4 ) 



v/8Asp^ 



(21) 



where 7(x) = y/irxe^ (1 — Erf (a:)) , Erf{x) is the error 
function, and 7 < 1 in this limit. By Eq.lj^OJ), 




FIG. 2: The geometric phase of the spin- 2 particle versus sjv 
and A. The spin was dissipatively coupled to an environment 
and the phase was calculated at time r = 27r/A, depending 
on the time scale of free evolution of the particle. 



it is clear that $g = in this limit, too. The numerical 
result for the geometric phase as a function of A and sm 
was shown in figure 2. To plot this figure, we assume 
that the system has evolved for time r = 27r/A, which is 
the characterized time for the system undergoing a free 
evolution. Figure 2 shows that the geometric phase is 
zero in the two limiting cases A 3> sat and A ^ s^r 
as expected. There are sharp changes among the line 
A ^ Sat, indicating a crossover from the limit A ^ sa? 
to A 3> Sat. 



V. SUMMARY AND DISCUSSION 

We have presented a general calculation the geometric 
phase in open systems subject to dephasing and dissipa- 
tion, the calculations are beyond the quantum Markov 
approximation and the weak coupling limit. For the 
dephasing system, it acquires no geometric phase with 
the decoherence rate 7 — > c», this can be explained as 
an effect of decoherence on the geometric phase, i.e., 
the quantum system could not maintain its phase in- 
formation after the decoherence time. There is a sharp 
change along the line 9 = 7r/2 as figure 1 shown, this 
can be understood in terms of the Bloch sphere that rep- 
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resents the state of the system. The geometric phase 
increases due to decoherence when initial states fall onto 
the upper semi-sphere, but it decrease when the ini- 
tial states on the lower semi-sphere. These results are 
similar to the prediction given by the quantum trajec- 
tory analysis for dissipative systems. The geometric 
phase $' in dissipative systems is always zero as long 
as Px{t)/py{t) =constant, this is exactly the case when 
A/sn — *■ oo or A/sn — *■ 0. A/s^ — > oo implies that 
the self energy A of the particle is much larger than the 
cumulative variance spf of the coupling constants ■ For 
gk taking the value +g or —g (g arbitrary) with equal 
probability, = J2k 9k - This tells us that the geometric 
phase is zero when the self-Hamiltonian dominates. On 
the other hand, when A/sn — > the interaction Hamil- 
tonian dominates, pointer states in this situation coincide 
with the eigenstates of the interaction Hamiltonian, thus 
the spin-i particle could not acquire geometric phase. In 
the crossover regime A ~ s^r, the geometric phase change 
sharply due to the interplay between the self-Hamiltonian 
and the interaction with the environment. 

These results constitute the basis of a framework to 



analyze errors in the holonomic quantum computation, 
where two kinds of errors are believed to affect its per- 
formance. This first error would take the system out of 
the degenerate computation subspace, while the second 
takes place within the subspace. The first kind of error 
can be eliminated by working in the ground states and 
having a system where the energy gap with the first ex- 
cited state is very large. The second kind of error falls 
to the regime analyzed in Sec. II and III, since there 
is no dissipation but dephasing in the system, while the 
first belongs to the regime discussed in Sec. IV. The cal- 
culation presented here in principle allows one to study 
the geometric phase at any timescale, and hence it has 
advantages with respect to any treatment with approxi- 
mations in most literatures. 
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